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The nonlinear c-model is considered to be useful in describing hadrons 

(Skyrmions) in low energy hadron physics and the approximate behavior of the global 

texture. Here we investigate the properties of the static solution of the nonlinear a- 

model equation coupled with gravity. As in the case where gravity is ignored, there 

is still no scale parameter that determines the size of the static solution and the 
re- 
winding number of the solution is 1/2. The geometry of the spatial hyperspace in 

<^ ' 

the asymptotic region of large r is explicitly shown to be that of a flat space with 

O ' some missing solid angle. 
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When a certain symmetry of a system is broken spontaneously, a masseless scalar particle 
naturally emerges as a Goldstone boson [|I[ for each broken generator of the symmetry The 
simplest lagrangian density, which has been used to demonstrate the above formal argument, 
can be written as 

c = U^ a d,r-K\^\ 2 -^? (i) 

where a = 0, 1,2,3 are chosen as an example in this paper. The physical vacuum of the 
system can be defined as a field configuration which satisfies | | 2 = v 2 . If we choose the 
0°-direction as the one that acquires the non-vanishing vacuum expectation value such that 

(<P°) = v 

(0*>=O, 2 = 1,2,3 (2) 

then the symmetry of the system, 50(4), is reduced to 50(3) and J, s become masseless 
Goldstone bosons. In the limit, A — > oo, the massive mode in the 0°-direction becomes 
infinitely heavy and decoupled from the system. The symmetry now can be realized in 
terms of J, s only: the unbroken symmetry is realized linearly but the transformation in the 
direction of broken generators nonlinearly. 

In terms of Goldstone bosons, the simplest lagrangian density in which the whole sym- 
metry of the system is properly incorporated as discussed above is the nonlinear a-model 
lagrangian density 

f 2 
£ = _^ Tr (^ f/ ^ f/ t) (3) 

where 

U = exp(ir-f/U) (4) 

The classical soliton of Eq. (|3|) has been found to be very important because of its use- 
fulness in describing hadrons (Skyrmions) 0] in low energy hadron physics as well as the 
approxmate behavior of the global texture which is attracting interest as a potential source 



of density perturbation for large scale structure of the universe || . Since the simplest form 
of Eq.(|3|) cannot support a stable soliton with finite energy ||, higher derivative terms are 
needed in Eq.([I|) or Eq.([3]) to describe the system with finite energy. They are the Skyrme 
terms used in low energy hadron physics. 

It is observed || that there is a set of solutions which extremize the action with the 
simplest lagrangian density, Eq.(f|), using the hedgehog ansatz: 

U = exp(zr ■ fF(r)) (5) 

The profile function F(r) satisfies the following equation; 

d 2 F 2dF 1 

+ z^n = -2^ 2F (6) 



dr 2 r dr r 2 

In this equation, one can see that there is no scale parameter which can determine the size 
of the soliton. This leads to the conclusion that, due to Derrick's theorem, there is no stable 
solution with finite energy. Fig.l shows a numerical solution of Eq.(^) with the boundary 
condition 

F(0) = vr. (7) 

It was analytically proved in Ref. || that there are no other solutions except for those 
solutions that can be obtained through the following transformations; 

F{r) — ► F{\r) (A = constant) 

— > F(r) + mr (n — integer) 
- -F(r) (8) 

The profile function F(r) in Fig.l approaches tt/2 as r — > oo as can be seen in the asymptotic 
form: 

^M-f + v^cosf^lnft+a) (9) 

2 V r \ 2 r J 

The winding number defined by 



= F(0)-F(oo) _ sin 2F(0) - sin 2F(oo) 

7T 27T l ' 

can be calculated to be \. Although the energy of this solution is infinite, it cannot be 
excluded mathematically as one of the solutions that extremize the action. It might be useful 
in describing the physical situations where a very large number of particles are involved in 
heavy ion collisions or where the solution has to be truncated due to the finite size of the 
correlated volume which is restricted by the horizon in the early universe. In Ref. 0, the time 
evolution of the textures whose initial configurations are given by small perturbations from 
the solution in Fig.l was investigated. Within the one parameter set of initial configurations 
considered there, it was found that textures with winding number larger(smaller) than 1/2 
collapse (expand). The critical nature of the winding number 1/2 was also discussed in 
Ref. J7|. The author presented the analytical predictions about the effect of the winding 
number on behaviors(time evolutions) of the scalar field and showed configurations of the 
predictions by numerically solving the field equations with a couple of particular forms of 
initial configurations. Although the above findings do not provide the full general proof, they 
lead one to suspect that the critical property that determines the collapse is whether the 
winding number is larger than 1/2. Ref. |8j presents numerical investigations of the collapse 
and unwinding of global textures without recourse to the nonlinear a-model approximation. 
As was explicitly shown in Ref. ||, the static solution shown in Fig.l is not a stable one. 
However, being a unique static solution that extremizes the action, it may represent the most 
probable initial configuration of texture formation. A small perturbation will eventually 
cause it to collapse or expand, but its collapse or expansion will be delayed longer than that 
of any other initial configuration. 

When we consider the situation with very large number of particles or with very large 
amount of energy, the effects of gravity and/or quantum corrections which can be put into 
the lagrangian in the form of higher derivative terms may become important. Hence it is 
interesting to investigate those effects on the winding number one-half solution. Here we 
consider the effects of gravity starting with the following action; 



/ = J d'x^Tg (t^R ~ f -fTr(d»Ud^)\ (11) 

with U given by Eq.(|J). For earlier works concerning the system of texture coupled with 
gravity, see for example Ref.'s 0, |TD| and [|TT|1 . 

For a static and spherically symmetric solution, the metric can be written as 

ds 2 = e 2q{r) dt 2 - e 2p{r) dr 2 - r 2 (d6 2 + sin 2 Od^). (12) 

The matter field equation obtained from 51 /5F = is 

d 2 dF dq dp 2 dF 2 . 

-r( r -T-) + \-r ~ -r) r ~r ~ e sin(2F) = (13) 

dr dr dr dr dr 

The equations for p(r) and q(r) obtained from the Einstein's equations, G^ v = —8ttG T m ", 
are 

where k = 8nGf 2 . From Eq's.([L3|), (p!4|) and ( [15| ), we can see there is no scale parameter, 
and hence no scale to determine the size of the solution. It is then an interesting question 
whether there exists a static solution which extremizes the action, Eg. (|TTD , which corre- 
sponds to Eq.(^) without gravity. In our calculations, we first put Eq.([l|) into Eq. (|1"3|) to 
eliminate q(r), couple the resulting equation with Eq.(|14D and numerically integrate them 
to obtain F(r) and p(r) with the boundary condition 

F(0) = Ti, p(0) = 0. (16) 

The only other input in the calculation is the value of dF/dr at r = 0. Changing the value 
induces only the change of the overall scale of the solution. The result of the calculation is 
given in Fig.2. Here also a series of solutions can be generated from the solution in Fig. 2 
through the transformations similar to Eq.(|8|), 



F(r),p(r) — »■ F(Xr),p(Xr) (A = constant) 
— >■ -F(r) + nir,p(r) (n = integer) 
--F(r),p(r) (17) 

Our numerical analysis indicates that 

| F(0) - F(oo) |= |, (18) 

although we still lack the analytic proof. The asymptotic forms at large r can be calculated 
explicitly to be 

F ^ o + V C0S o ln — + 7 ' 19 

2 V r yi r y 

p -> o + 6— sin (/31n(-) + 7 + s\ , (20) 



g — >• c - (1 + — - J- ———cos pln{—) + 27 



/3 d kx 



4(1 + /3 2 ) 4r 



sin (/3ln(—) + 2 7 ") - b-sin (pin(—) + 7 + s) (21) 



where 



r - ^. (23) 

2 = _+_^ (24) 

2561 -«' V ; 

tan 5 = -P, (25) 

and r , 7 and c are integration constants. We can put c to be zero by adjusting the scale of 
the coordinate t so that it coincides with the scale of proper time in the asymptotic region 
of r — ► 00. 

Let us now consider the geometry of the hypersurface of t = constant in the asymptotic 
region of large r. The asymptotic form of the metric is given by 

ds 2 = -^—dr 2 + r 2 (d9 2 + sin 2 Odtf) (26) 

1 — K 



or with the coordinate transformation, r = Jj^r, 

ds 2 = df 2 + (1 - K)r 2 {dd 2 + sin 2 6d(j) 2 ) (27) 

This metric represents a flat space in the asymptotic region with the deficit solid angle 

5Q = 47tk (28) 

To see that the metric represents a flat space in the asymptotic region, we calculate all 
components of the Riemann curvature tensor in the orthonormal basis which we take as 



e 1 = J- dr, e 2 = rd6, e 3 = rsin9d(f) (29) 



The only nonzero components are 



R 2 ™ = -£ (so) 



and those related to it by algebraic symmetry or antisymmetry. We see that all components 
of the Riemann curvature tensor in the orthonormal basis converge to zero as r — > oo. It 



has been shown in Ref. [O that the self-similar texture solution also produces a flat spatial 
hypersurface with some missing solid angle in the asymptotic region of large r. 

We next consider the energetics of our solution. With the use of the hedgehog ansatz 
(Eq.(||)) and the metric form of Eq . (|12"D , the tt-component of the energy-momentum tensor 
T^ v is calculated to be 

T "4f^ 2! '(^) ! + e "^}- < 31 > 

To investigate the behavior(r - dependence) of the energy density in the asymptotic flat 



region of r — > oo, we put the asymptotic solutions, Eq.'s (|19D, ([20D and ([2lD , into Eq. (|3lD 
and obtain the result 

T u ~ — (32) 

as r — ► oo. This shows that the total energy contained in a volume diverges linearly in r as 
the volume is increased. 



In conclusion, we see that the basic structure of the classical static solution is not changed 
by gravity: there is still no scale parameter that determines the size of the field configuration 
and the solution represents the knot of windung number 1/2. The geometry of the spatial 
hypersurface in the asymptotic region of large r is that of a flat space with some missing 
solid angle. 
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FIGURES 
FIG. 1. The profile function F(r) without gravity. 

FIG. 2. The profile function F(r) and the function p(r) with gravity (« = 0.865). 
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http://arXiv.org/ps/hep-th/9407072vl 



s 



2 - 
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